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Abstract 

We map the complex electric fields associated with the scattering of surface plasmon polaritons 
by single sub- wavelength holes of different sizes in thick gold films. We identify and quantify 
the different modes associated with this event, including a radial surface wave with an angularly 
isotropic amplitude. This wave is shown to arise from the out-of-plane electric dipole induced in 
the hole, and we quantify the corresponding polarizability, which is in excellent agreement with 
electromagnetic theory. Time-resolved measurements reveal a time-delay of 38 it 18 fs between the 
surface plasmon polariton and the radial wave, which we attribute to the interaction with a broad 
hole resonance. 

PACS numbers: 73.20.Mf, 78.68.+m, 73.50.Bk 
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The interaction of light, and in particular of surface waves such as surface plasmon po- 
laritons (SPPs), with sub- wavelength holes in metallic films leads to phenomena such as 
extraordinary optical transmission (EOT)^'^ and negative index metamaterials^^"-. In this 
context, Sub- wavelength holes are important in diverse fields of research, such as nanopho- 
tonics, enhanced nonlinear optics, and biosensing, and consequently, with a few notable 
exceptionsiiii^, experiments focus on the macroscopic consequences of the light-hole inter- 
actions. Conversely, much theoretical work has been devoted to their understanding from 
a microscopic viewpoint^"—. That is, while the models intrinsically consider scattering of 
light and SSPs from sub-wavelength defects, often in terms of their polarizabilities, most 
experimental studies address collective response features such as EOT spectra and their 
dependence on changes in sample geometry, film thickness^, and hole shaped. In essence, 
a detailed investigation of the dynamics of plasmonic scattering remains incomplete. 

In this Letter, we present a systematic experimental study of the scattering of SPPs from 
individual circular sub-wavelength holes of different diameters in optically thick gold films. 
We use phase-sensitive time-resolved near-field optical observations to measure the plasmonic 
scattering on a nanometer length- and femtosecond time-scale, determining the angularly 
resolved scattering amplitude for sub- wavelength holes of different sizes and quantifying the 
temporal dynamics of the scattering events. We show that these results can be explained in 
terms of the polarizability of the hole, in excellent agreement with electromagnetic theory. 

We investigate the scattering from individual holes, with diameters ranging from 50 to 
980 nm, milled with a focused ion beam into a 200 nm thick gold film deposited on a 
glass substrate. We launch SPPs on the air-gold interface by illuminating a slit in the film 
with either 150 femtosecond FWHM pulses or a continuous wave laser source centered at a 
wavelength of 1550 nm. At this wavelength, a waveguide mode is supported by holes with 
diameters above 760 nm. These SPPs are efficiently directed towards the hole by a Bragg 
grating that is milled into the film on the opposite side of the slitr^, as shown in Fig. [H 
We measure the complex electric field of this system, at a height of 20 nm, with near field 
microscopy^. 

In Fig. |2] we present typical experimental results for a hole with a diameter of 875 nm. 
We observe fringes in the amplitude that originate from the interference of the forward 
propagating SPP with the scattered waves. Here, the dark vertical stripes near 5 and 8 
fim are attributed to geometrical imperfections in the area of the launching slit. The salient 
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FIG. 1: (Color online) Scanning electron micrograph of a SPP launcher and a hole. The arrow 
shows the direction of propagation of the excited SPPs. Illumination is provided from underneath 
the film. 

features of the data are reproduced by first-principle calculations [Fig. Mjo)] performed using 
a modified version of a model presented elsewhere^^. This method, which involves a rigorous 
expansion of Maxwell's equations, both in terms of plane waves on either side of the film 
and in cylindrical waves inside the hole, is used to calculate the electric and magnetic 
polarizabilities of the hole, which are in turn used to obtain the scattered fields^. 

The different waves present in the scattering event are readily identifiable in k-space^ 
after a Fourier transform [Fig. [2tc)]. In A;-space, signatures of the forward propagating SPP, 
evidenced by the peak near 2°, and the scattered wave, which emerges as a circular pattern, 
are readily separable. These signatures are shown in Fig. Et^d), in which both the data (thick 
shaded curves) and the fits for the different surface waves (thin dark curves) are presented. 
The distributions of both waves peak around a wave vector k = 4.2 /im~^, which is in good 
agreement with the theoretical value of kgp = k^^e/ (e + 1))°'^ = 4.1 /xm"^, obtained from 
the dielectric constant of gold at 1550 nm, e = — 115 + 11*^^. 

In the inset of Fig. [2|^d), we show the angular dependence of the in-plane electric field 
distribution. Here, the outer curve represents the total integrated field amplitude while the 
inner curve represents the integrated amplitude that we can attribute to the surface waves. 
We find these curves by integrating up to /c = 5.0 /im"^ under the data and fits, respectively, 
in each direction (main part of the figure). 

First, we observe no angular dependence of the radial mode and attribute the slight left- 
right asymmetry found in our measurements to either a slightly asymmetric hole or near-field 
tip. Second, for this relatively large hole the amplitude content of the radial wave is about 
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FIG. 2: (Color online) Typical near- field measurements of plasmon scattering from a sub- 
wavelength hole, (a) Real space image at the time when the peak of the SPP pulse is in the 
scan area. The pulse has gaussian profile and is incident from the bottom of the image. The hole 
position is indicated by a solid circle. Fringes caused by interference between incident and scattered 
plasmons are clearly visible, (b) Calculated electric field amplitude under the conditions of (a), 
(c) A;-space representation of the data shown in (a). Incident and transmitted SPPs are associated 
to the vertical line, while scattered plasmons appear as a circle. The image has been saturated to 
make the signature of scattered modes clearly visible, (d) Cuts along 2° and 35° radial directions 
in A;-space [see arrows in (c)] showing signatures of forward and radial plasmon waves, respectively. 
The 35° line has been magnified by a factor of 70, its amplitude offset by 5, and position shifted 
by 2 /im~^ for clarity. Thick and thin curves are the data and the corresponding fits, respectively. 
Inset: angle-dependent integrated field amplitude; the outer blue curve is integrated up to /c = 5 
fim~^ (see dashed curves in main plot); the inner purple curve represents the combined amplitude 
of the two surface modes identified in (c). 
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7% of the amplitude of the forward propagating SPP. As we show below (Fig. |5]), for an 
875 nm hole ~ 0.7 of the forward propagating SPP amplitude is contained in the incident 
beam, and hence the amplitude of the radial wave is ~ 10% that of the incident SPP. Note 
that our incident Gaussian beam is ~ 6 fim width, which is larger than the hole (see^^ for 
a discussion on the scattering cross section of the hole. Finally, from this data we can also 
estimate that approximately half of the radially scattered field amplitude is in the surface 
modes. The remaining half of the radial fields have in-plane wave vectors that are smaller 
than ko = 27r/A = 4.05 yum"^, corresponding to radiative modes of light, which propagate 
upward, away from the metal surface. Qualitatively, these are seen as the regions of non-zero 
amplitude within the surface wave circle in Fig. [2]^c), or in the difference between the inner 
and outer curves (other than near 0°) in Fig. Wid). 

In order to determine the nature of the observed scattered field, we consider its asymptotic 
formic 

X e-if e'('»«+"'»=)/l" {^) (-uvfl + , (1) 

where the azimuthal scattering amplitude 

A*"^ (ip) = cos (p — vT+e ruy cos (p + ^/e (2) 

is given for excitation by a SPP that travels along the x direction. Here, R is the in-plane 
radial distance to the hole center, z is the height above the film surface, Wgp = —ko/ (e + l)'^'^ 
is the out-of-plane component of the SPP wave vector, z and R are out-of-plane and in-plane 
unit vectors, px and Pz are the in- and out-of plane electric dipoles, and my is the in-plane 
magnetic dipole of the hole. Consequently, the isotropic nature of the detected scattered 
field [inset to Fig. EJj^d)] suggests that we are mostly sensitive to the contribution from p^. 
This is not entirely unexpected since the contributions of other dipole components peak near 
ip = 0° and 180° and are therefore hard to separate from the incident SPP field. However, 
we note that to accurately reproduce the field in all directions [Fig. MJo)] the fields that arise 
due to well as those due to Pz, are required^. 

We image the pz dependent electric field by filtering out the forward propagating SPP 
mode in /c-space and then Fourier transforming back to real-space. In Fig. [3] we show this 
field, both near cp = —90° where only a contribution from pz is expected, and in a 20 /xm^ 
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FIG. 3: (Color online) Radial distribution of the field scattered by a 875 nm hole. The measured 
real part of the electric field is shown by dots, while the curve is an exp [ikR)R~'^ fit. The inset 
shows the two-dimensional image of the real part of the electric field, and the arrow indicates the 
line along which the data shown in the main part of the figure is taken. 

area around the hole (inset). A fit of the form exp {ikR)R~°' reproduces the experimental 
data with k = 4.20 ± 0.02 fim~^ and a = 0.57 ± 0.10. Hence, the measured amplitude 
damping is in excellent agreement with the l/y/R behavior predicted for circular waves 
by Eq. ([1]), as required for energy conservation for radial propagation on a plane. The 
discrepancy between the data and fit close to the hole may be indicative of the presence of 
creeping waves, which can have substantial amplitudes at short distances^. We can also use 
Eq. dl]) to express the ratio of the scattered field to the incident field for ip = 90° in terms 
of the electric polarizability 



Hence, we measure = (4.1 ± 0.1) x 10~^^ m^, which is larger than, but comparable to, 
the value |q!£;| = 3.3 x 10~^^ m^ calculated for a hole in a perfectly conducting filn>i^. 

In Fig. m we present the temporal dynamics of plasmonic scattering, showing the ampli- 
tudes of both the forward and scattered radial modes as a function of time. The forward 
amplitude is averaged over a range of ±5° while the scattered amplitude is averaged over 
a similar range in the transverse directions {6 = ±90°) to ensure no overlap between the 
modes. Both amplitudes are normalized to the peak of the forward propagating SPP, and 
the second peak of the double peak results from SPPs that are initially launched backward 
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FIG. 4: (Color online) Temporal dynamics of plasmon-hole scattering. The normalized forward 
field amplitude (circles, left axis) and the corresponding fractional scattered field (squares and 
diamonds, right axis) are shown as functions of the time delay. The solid curves are double-peak 
fits which provide guides to the eye; the dashed curves show the individual peaks of the scattered 
amplitude fit; these lag behind the peaks of the forward amplitude (vertical dotted line at 700 fs, 
for the latter peak) by ~ 40 fs. 

and subsequently reflect from the Bragg grating towards the hole^^. Most strikingly, we 
observe that the radial wave peaks one time-step after the forward propagating SPP: there 
is a time-delay between the forward propagating SPP pulse and the radial pulse resulting 
from the plasmonic scattering. By averaging over the time-delay that we measure between 
incident SPP and the scattered wave, in all available directions, we quantify this shift as 
38 ± 18 fs. While the frequency dependence of the polarizability of the hole alone results in 
a shift only a few fs, it is able to explain the observed shift when combined with a modestly 
chirped incident pulse. That is, since the different plane waves that form our pulse scatter 
from the hole with different amplitudes and phases, their reconstitution into a beam results 
in the observed time shift. This shift, then, represents a first measurement of the complex 
spectral response of an individual hole. 

We study the effect of hole size on the plasmonic scattering for hole sizes ranging from 50 
to 980 nm. The results are summarized in Fig. [5l The inset depicts the fraction of the field 
amplitude that is contained in the forward propagating SPP mode, which we normalize to 
the field on films without a hole. As expected, less energy is found in this plasmonic wave 
as the hole size increases and hence scatters more energy. 
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FIG. 5: (Color online) Hole-size dependence of the plasmonic scattering. We plot the radially 
scattered field amplitude normalized to the peak of the forward propagating SPP as a function 
of hole size. The solid curve shows the calculated hole-size dependence of pz (and hence the 
calculated scattered amplitude). Inset: fraction of field contained in the forward propagating SPP 
peak, normalized to unity for a featureless film, as a function of hole size. The vertical dashed lines 
show the cutoff of the lowest-order guided mode in a deep hole for 1550 nm light. 

The main part of Fig. depicts the dependence of the scattered field amplitude of the 
radial wave normalized to that of the forward plasmon field (symbols) on the hole diameter, 
d. This dependence is accurately reproduced by the calculated electric polarizability of the 
hole (curve )^. An offset of 0.017 is present in the data (and added in the calculations) due 
to measurement noise, and represents a signal-to-noise ratio of about 50 to 1, while the error 
bar is mainly indicative of the asymmetry of our measurements. For the smaller holes, well 
below the cutoff, the field amplitude increases as ~ d'^, where the divergence from a perfect 
dipolar behavior is due to the finite size of the hole and the film thickness. As the hole 
size increases, and in particular once the cutoff is reached and waveguide modes in the hole 
become accessible, the scattered amplitude begins to saturate. We attribute this behavior to 
an increase in energy flow through the hole^^. Finally, we note that the scattered amplitude 
peaks at about 10% of the forward propagating SPP amplitude. 

To summarize, we have captured the near-field interplay between the waves associated 
with plasmon scattering from sub-wavelength holes and mapped both their dependence on 
hole size and their temporal dynamics. We show that both the hole size dependent amplitude 
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and the time-delay of the scattered mode can be understood in terms of the calculated po- 
larizability of the hole. Our results provide a comprehensive understanding of hole-plasmon 
scattering and should, in the future, allow for the optimization of the interaction between 
holes. Consequently, the geometric properties of hole arrangements can be tuned to cre- 
ate electromagnetic hot spots. These can be used to, for example, enhance emission from 
dye molecules, address individual quantum emitters, or create near field distributions for 
sensing or imaging purposes. Hence, these results have implications for a broad range of 
nanoplasmonic applications. 
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